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Abstract 

The Earth’s Moon is thought to have formed by an impact between the Earth and an impactor around 4.5 billion 
years ago. This impact could have been so energetic that it could have mixed and homogenized the Earth’s mantle. 
However, this view appears to be inconsistent with geochemical studies that suggest that the Earth’s mantle was 
not mixed by the impact. Another plausible outcome is that this energetic impact melted the whole mantle, but the 
extent of mantle melting is not well understood even though it must have had a significant effect on the subsequent 
evolution of the Earth’s interior and atmosphere. To understand the initial state of the Earth’s mantle, we perform 
giant impact simulations using smoothed particle hydrodynamics (SPH) for three different models: (a) standard: a 
Mars-sized impactor hits the proto-Earth, (b) fast-spinning Earth: a small impactor hits a rapidly rotating proto-Earth, 
and (c) sub-Earths: two half Earth-sized planets collide. We use two types of equations of state (MgSiOs liquid and 
forsterite) to describe the Earth’s mantle. We find that the mantle remains unmixed in (a), but it may be mixed in (b) 
and (c). The extent of mixing is most extensive in (c). Therefore, (a) is most consistent and (c) may be least consistent 
with the preservation of the mantle heterogeneity, while (b) may fall between. We determine that the Earth’s mantle 
becomes mostly molten by the impact in all of the models. The choice of the equations of state does not affect these 
outcomes. Additionally, our results indicate that entropy gains of the mantle materials by a giant impact cannot be 
predicted well by the Rankine-Hugoniot equations. Moreover, we show that the mantle can remain unmixed on a 
Moon-forming timescale if it does not become mixed by the impact. 
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1. Introduction 


The so-called giant impact hypothesis is a widely ac¬ 


cepted explanation for the origin of the Moon (|Hart- 


|mann and Davis|p^975| [Cameron and Ward[|197^ . Ac¬ 

cording to the standard version of this hypothesis, a 
Mars-sized impactor hit the proto-Earth and created a 
disk around the planet from which the Moon accreted. 
This hypothesis has been favored because it can ex¬ 
plain the Moon’s mass, iron depletion, and the angular 
momentum of the Earth-Moon system. However, this 
model has difficulty in explaining the fact that the Earth 
and Moon have nearly identical isotopic ratios (e.g. 
oxygen, silicon, and tungsten, Wiechert et al.|200l| Her- 


wartz et al.|2014t[Armytage et al.|2012t|Touboul et al. 


20071. The typical outcome of a giant impact simula¬ 
tion is that the disk materials are derived mainly from 


the impactor (e.g., Benz et ak) 1986| ), which is often as¬ 
sumed to have had different isotopic ratios given that 
the oxygen isotopic ratios between the Earth and Mars 
differ by 0.321%o ( [Eranchi et aH |1999| ). [Reufer et al.| 
( 2012| ) report that an impact at a high impact velocity 
and steep impact angle would deliver more Earth’s man¬ 
tle materials to the disk, but it is still difficult to explain 
the identical isotopic ratios. 

Pahlevan and Stevenson ( |2007| ) have suggested that 


turbulent mixing in the Earth’s atmosphere and the disk 
homogenized the isotopic ratios of the two reservoirs. 
This model could potentially solve the isotopic problem, 
but it has several shortcomings. This mechanism may 
not work for all of the observed isotope systems, such as 


Si ( [Pahlevan et al.[[2011t|Armytage etn^|2012| ). Eur- 
thermore, this mixing would have required the Earth’s 
whole mantle convection during the Moon formation, 
but such convection may not have occurred due to the 
thermally stratified structure of the mantle after the im- 
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pact (discussed in Sections 2.4 and 3.3). Even if the 
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post-impact mixing caused the disk to have the same 
isotopic reservoir as, say, the outer 80% of the man¬ 
tle but failed to equilibrate with the inner 20%, then 
the Earth and Moon could still be isotopically differ¬ 
ent if there had been a subsequent mixing of the Earth’s 
mantle after the Moon formation. At present, there 
is no detectable oxygen isotopic difference (with re¬ 
spect to the three isotopes, which are referred to as 
- 0.52^^^O) among Earth rocks. 

The possibility that the impactor had a similar oxy¬ 
gen isotopic ratio to the Earth has been recently re¬ 
visited. |Kaib and Cowan| ( |2015| ) investigate feeding 
zones of terrestrial planets based on orbital calcula¬ 
tions and predict that this possibility is ~ 5% or less. 
Mastrobuono-Battisti et al.| ( [20151 ) have performed sim¬ 


ilar analyses and report that the possibility could be as 
high as 20 - 40% if its standard deviation (±cr) is in¬ 
cluded, although the most plausible value of this possi¬ 
bility (~ 10 - 20%) may not be significantly different 
from the value estimated by |Kaib and Cowan] ( [2015| ). It 
is also possible that the impactor could have been com- 


positionally similar to enstatite chondrites (e.g., Her- 


|wartz et SI1|2014 ), whose compositions are much more 
similar to those of the Earth than those of Mars. Al¬ 
ternatively, a recent model of planet formation, the so- 
called Grand Tack model ( Walsh et^|2011 1, may sug¬ 
gest a different outcome. This Grand Tack model sug¬ 
gests that the planetesimal disk was truncated at 1 AU 
due to migration of gas giants and, as a result, terrestrial 
planets mainly formed from the inner part of the disk. 
This model may increase the chance of having an im¬ 
pactor with a composition similar to the Earth (personal 
communications with Alessandro Morbidelli). This in¬ 
creased probability may arise possibly because the main 
source of the Earth’s materials was confined to a limited 
region of distances from the Sun, or perhaps because of 
more efficient mixing of the source materials than the 
standard models predict (discussed in Section 4.4). In 
either case, this finding would imply that the Earth is 
different from Mars but not necessarily different from 
the terminal giant impacting body that led to the forma¬ 
tion of the Moon (often called as “Theia”). It should be 
noted that having the same oxygen isotopic ratios for the 
proto-Earth and impactor does not necessarily explain 
the nearly identical tungsten and silicon isotopic ratios 


of the two (Dauphas et al. 2014 Eitoussi and Bourdon 

|20T2l ). 

New gian t impact models have bee n suggested as al 
ternatives. 


Cuk and Stewart 


(2012) propose that an 


impactor hit a rapidly rotating proto-Earth (called the 
“fast-spinning Earth”), whereas Canup ( 2012| ) suggests 
a giant impact between two half-Earth-mass objects 


(here we call this “sub-Earths”, and the model is also 
called “large impactor” and “half Earth” in other lit¬ 
eratures). In these cases, the composition of the disk 
would have been similar to that of the Earth’s mantle; 
therefore, the models could explain the isotopic simi¬ 
larities. In these models, the angular momentum of the 
Earth-Moon system a fter the impact was 2-3 times as 
large as today’s value. 


Cuk and Stewart 


(2012) suggest 


that the evection resonance between the Moon and the 
Sun could have transferred the excess angular momen¬ 
tum of the Earth-Moon system to the Sun-Earth system. 
This resonance occurs when the period of precession of 
the pericenter of the Moon is equal to the Earth’s or¬ 


bital period ( |Touma and Wisdoni||1998| ). It is also pos¬ 
sible that there is some other resonance that yields the 
same result ( [Wisdom and Tian[ |2015j ). Nonetheless, it 
is not yet clear whether any resonance could efficiently 
remove the excess of the angular momentum (personal 
communications with Jack Wisdom). The existence of a 
resonance does not imply the removal of a large amount 
of angular momentum because that would depend on a 
particular and possibly narrow choice of tidal parame¬ 
ters. 

These new models are indistinguishable in terms of 
the oxygen isotopic ratios, but additional geochemical 
constraints may differentiate these models. Eor exam¬ 
ple, it has been suggested that the Earth’s mantle may 
not have been completely mixed by the giant impact. 
This has been drawn from various isotopic studies, es¬ 
pecially those on the Hf-W system. Hf is a lithophile 
(“rock-loving”) element, whereas W is a moderately 
siderophile (“iron-loving”) element. ^^^Hf decays to 
i^^w with a 9 Myr half-life; thus, the mantle of a planet 
would have an enhanced i^^W/^^'^W if differentiation 
occurred while Hf was still alive (within the first ~ 60 
Myr after CAI formation). Most terrestrial rocks hav e 
similar values of ( Lee and Halliday 1996| , 

but jWillbold et al.| ( [2011| ) and Touboul et al.| ( 2012| ) find 


that 2.8 and 3.8 billion years old rocks have ratios in ex¬ 
cess of these values. This finding may indicate that the 
early mantle was heterogeneous (while ^^^Hf was still 
present) and that the signature was preserved at least 
until 2.8 billion years ago. Determining the cause of the 
heterogeneity is an active area of research. It should be 
noted that the chemical heterogeneity could take many 
forms, including a non-uniformity of noble gas mole 
fraction (e.g., [Mukhopadhy^[2012] ) or a discrete layer 
of denser mantle material at the base of the mantle be¬ 
cause of the formation of a basal magma ocean (e.g., 
[Labrosse et al.|[2007] ). 

Regardless of the cause or form, if the mantle hetero¬ 
geneity formation predated the Moon-forming impact. 
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the Earth’s mantle may not have been mixed even by 
the giant impact. Previous studies based on the lunar 
I82w/i84w ratio suggest that the Moon formed as early 
as 30 Myr after CAI formation (e.g., |Lee et al.[pW7] ), 


whereas [Touboul efaL ( 2007| ) propose that ~60% crys¬ 
tallization of the lunar magma ocean occurred after ~60 
Myr by taking into account the excess of formed 
by neutron capture of ^^^Ta. This is consistent with the 
age estimate based on other systems (e.g., Sm-Nd, [Carl-| 
son and Lugmair|1988 and recent studies on orbital dy- 


Jacobson et al.|2Q14| ). Although the age of the 


namics, 

Moon is still under debate, recent work tends to suggest 
a young age of the Moon (> 60 Myr). Herein, we fo¬ 
cus on the scenario in which the mantle heterogeneity 
predated the formation of the Moon. 

In addition to the mantle mixing, understanding the 
extent of the impact-induced mantle melting is impor¬ 
tant because the extent affected the evolution of the 


Earth’s interior and atmosphere (e.g., Abe and Matsui 
|1986t[Tbnks and Melosh[ p^93l |Elkins-Tanton| 2008| ). 

When the Earth grew through collisions with numerous 
impactors, these impactors melted part of the Earth’s 
mantle and delivered their metallic iron to the Earth. 
The metallic iron of small impactors (at least up to a few 
hundred kilometers in size) would have been dispersed 
as droplets and the resulting rainfall would have led 
to metal-silicate equilibration during the descent to the 
metal pond at the base of the magma ocean ( [Stevenson 


1990). The iron might have passed through the solid- 


rich deeper mantle without further equilibration with 
the ambient mantle. If the abundance of siderophile el¬ 
ements in the mantle recorded the metal-silicate equi¬ 
librium at the base of the magma ocean, the magma 
ocean depth would have been approximately around 28- 


40 GPa (700-l,200km deep) (e.g., Li and Agee 1996 


[Righter et al.j |1997| ). However, this model may be too 


simplistic because it assumes that the core formation oc¬ 
curred by a single stage process, but the Earth’s core 
must have formed though multiple impacts processes 
( jWade and [200^ [RUbie et al.j [2^ . Thus, the 

mantle geochemistry is not likely to have recorded the 
single impact event. The concentrations of siderophile 
elements in the mantle reflect very complex processes 


man et al.||1999l|Rubie et al.||2003l|Dahl and Steven- 

sonl|2010l|Shietal.l|2013|l. 


Analytical and numerical studies suggest 
nificant fraction of the mantle would have e 
melting by the Moon-forming impact (e.g.. 

that a sig- 
xperienced 
Tonks andj 

Melosh 1993[ Canup||2008[ de Vries et al. 

2014). A 


simple estimate can be described as follows; for an im¬ 
pacting body with velocity 10 km/s, the specific kinetic 


energy carried by the body is 5 x 10^ J/kg. The la¬ 
tent heat of melting is about 1 x 10^ J/kg (for a mean 
temperature of 2500K). Therefore, a Mars-mass projec¬ 
tile would deliver several times more energy than that 
needed to melt the entire mantle, assuming that the pre¬ 
impact state is near the solidus. Of course, this does not 
indicate that the entire mantle is in fact melted upon im¬ 
pact because the heating is heterogeneous and because 
part of the impact-induced energy is partitioned to the 
kinetic energy of the system. This makes the impor¬ 
tant point that the extent of melting must be assessed by 
considering the budgeting of all delivered energy and 
not merely by considering the shock heating associated 
with the primary impact (discussed in |4.1| ). In this re¬ 
spect, giant impacts differ from small impacts. 

It is often assumed that if there is complete melting 
at one time, then the preservation of geochemical reser¬ 
voirs or the retention of volatiles (e.g.. He) is not pos¬ 
sible. To the extent that this assumption is correct, the 
geochemical evidence appears to contradict models in 
which there is complete melting. However, we will ar¬ 
gue that preserving heterogeneity (i.e., lack of complete 
mixing) is possible even if most or all of the mantle be¬ 
comes molten by an impact (Section [2^. 

We also consider another problem: even if the man¬ 
tle heterogeneity survived the Moon-forming impact, 
it could have been erased in the subsequent evolution. 
This could have occurred in two different ways: in the 
period immediately after the giant impact during the pe¬ 
riod of “healing” (cool-down following a giant impact, 
leading to a thermal state rather similar to that before 
the giant impact); or during the very long evolution of 
the magma ocean and solidifying mantle (time scale of 
millions to even hundreds of millions of years). This 
paper is concerned only with the first of these. 

We aim to make connections between dynamics and 
geochemical observations in terms of the Earth’s early 
mantle. Herein, our main questions are as follows: (1) 
What is the extent of mantle melting after the Moon¬ 
forming impact? (2) If the pre-impact Earth’s mantle 
was heterogeneous, did the giant impact erase the sig¬ 
nature? (3) Was the post-giant impact cooling (during 
the subsequent thousands of years) capable of erasing 
the pre-giant impact heterogeneity? We attempt to an¬ 
swer these questions using the three different models: 
(a) standard, (b) fast-spinning Earth, and (c) sub-Earths. 
Smoothed particle hydrodynamics (SPH) is used for the 
giant impact simulations (Section |2.1|). Two different 
equations of state are used (Section p!^ . Eurthermore, 
we investigate the possibility that the Moon-forming 
impactor and Earth were isotopically similar (Section 
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2. Model 

2.7. Smoothed particle hydrodynamics (SPH) 

SPH is a Lagrangian method in which a fluid is rep¬ 
resented by numerous particles (grids). A particle i has 
a characteristic size hi, which is the so-called smooth¬ 
ing length. The mass of the particle m/ is distributed 
within a sphere of radius 2hi. Each particle has its own 
density distribution according to its weighting function 
(kernel). The density at a given location is calculated 
as the sum of the density distributions of nearby par¬ 
ticles. In a standard SPH, each particle has the same 
mass, therefore a massive object (e.g. a planet) can be 
resolved better than a less massive object (e. g., a Moon¬ 
forming disk). The details of SPH are included in pre¬ 
vious studies (e.g., |Nakajima and Stevenson||2014| ). 


2.2. Equation of state 

We use two different equations of state (EOS) for 
the mantle materials. One of the equations repre¬ 
sents forsterite Mg 2 Si 04 (hereafter “forsterite EOS”). 
Eorsterite is the magnesium end member of the olivine 
solid solution series and is the major mineral phase in 
the upper mantle (P<~14GPa). The forsterite EOS is 
derived from the semi-analytic equation of state known 
as M-ANEOS ( [Thompson and Lauson[ |1972t jMeloshj 
|2007| ). M-ANEOS has been commonly used for hy¬ 
drodynamic simulations and it can treat phase changes 
and co-existing multiple phases. However, it does not 
correctly describe the high-pressure properties, espe¬ 
cially when the starting material is already in a high- 
pressure phase such as perovskite MgSiOs (now more 
correctly designated bridgmanite). Eor this reason, we 
also use an equation of state that represents MgSiOs liq¬ 
uid (“MgSiOs liquid EOS”). Since the Moon-forming 
impact would induce global melting of the mantle as 
discussed in SectionjT^ the liquid EOS is more suitable 
for this calculation. Eor comparison, we perform a cal¬ 
culation using an MgSiOs bridgmanite EOS and discuss 
the results in Sections [T^ and [M3] (Supplementary ma¬ 
terials). Since there is no ready-to-use MgSiOs liquid 
EOS for hydrodynamic simulations, we reconstruct an 
equation of state based on an analytic expression of the 
Helmholtz free energy and thermal coefficients derived 


from the first-principles calculations (e.g., de Koker and 
Stixrudej |2QQ9| ). The details are described in Section 


STT] 


2.3. Melting criterion 

We consider that a region is molten if its entropy gain, 
which is the entropy difference before and after the im¬ 
pact, exceeds the entropy required for melting A^meit- A 


wide range of A^meit has been reported, but here we take 
a high end-member A^meit ~ 623 J/K/kg in the subse¬ 
quent calculations (MgSiOs, Stixrude and Karki|2005| ). 
We also assume that the mantle is close to the solidus 
before the impact, as discussed in detail in Section [TS] 
It should be noted that since the entropy of melting is a 
direct measure of the disorder accompanying melting, it 
should be relatively insensitive to pressure or tempera¬ 
ture, unlike the latent heat of melting. 

2.4. Mixing criterion 


As discussed in Section |3.2| in greater detail, one of 
the most important consequences of a giant impact is 
that the outermost part of the planet would undergo a 
greater rise in irreversible entropy than the deeper re¬ 
gions because it would undergo higher shock-heating 
on average. This leads to dS jdr > 0, where S is the en¬ 
tropy and r is the mean radial distance from the center of 
the Earth to an equipotential surface (the equipotential 
surfaces within the planet are oblate due to the fast rota¬ 
tion of the planet). Therefore, the mantle is stable with 
respect to thermal convection immediately after the gi¬ 
ant impact. Under this circumstance, even if the mantle 
is molten, the mantle would not spontaneously convect 
or homogenize. After a giant impact, there is a variation 
in the entropy rise within the same r, but we ignore the 


azimuthal dependence (discussed in Section 4.2). Our 
calculations are carried out using entropy values that 
are averaged over latitude and longitude. We discuss 
the possible shortcomings of this approach in Sections 
lOand lSTTl 

The key issue here is that this entropy profile 
(dS I dr > 0) does not necessarily mean that the man¬ 
tle remains unmixed. The giant impact transfers a 
large amount of kinetic energy to the mantle, which 
may dynamically mix this stably stratified mantle. Our 
SPH calculations dampen the velocity shears on a short 
timescale (but without allowing turbulence), indicating 
an unrealistically large viscosity, and therefore they are 
not suitable for analyzing the mixing directly. We must 
instead use the velocity fields and energy state at an 
early stage after the giant impact to compute the extent 
of mixing. 

We investigate the mixing state of the mantle based 
on an energy balance. We assume that the mantle is un¬ 
stable and mixed if the impact-induced kinetic energy of 
the mantle AKE exceeds the potential energy of the sta¬ 
bly stratified mantle APE by approximately a factor of 
two, i.e., APE/AKE < 0.5 (see Section [S1.2| ). The most 
important source of AKE is differential rotation arising 
from the oblique impact, which is a necessary attribute 
of any giant impact that is capable of forming the Moon. 
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AKE is calculated as the kinetic energy differences in 
the mantles between before and after the impact, which 
is directly obtained from the SPH calculation (here, “af¬ 
ter” indicates a state where the Earth rotates as a rigid 
body). The key point here is that AKE is not the same as 
the kinetic energy of the post-impact Earth if the Earth 
has an initial spin (this is the case for b). This is be¬ 
cause the velocity shear is induced by the impact, not 
by the initial rotation of the planet and it is this differ¬ 
ential rotation that is responsible for the instability. The 
angular momentum is of course conserved in this pro¬ 
cess of energy redistribution. In an overturn that leads 
to mixing, half of the AKE contributes to overcoming 
gravity, and half is dissipated as heat. Therefore, in this 
case, the AKE is sufficient, by more than a factor of two, 
to overcome the potential energy of the stably stratified 
mantle. 

APE is the potential energy difference between the 
stability stratified mantle {dS/dr > 0) and mixed man¬ 
tle (S(r) = Aave =const., where Aave is the averaged en¬ 
tropy of the pre-mixing mantle). APE can be expressed 
approximately as 


4nr^V(r)Apdr, (1) 

^CMB 

where 

Ap = (^) (5(r)-5ave)- (2) 

^CMB is the distance from the center to the core-mantle 
boundary (CMB) and is the planetary radius. V(r) is 
the potential energy per unit mass at r. The gravitational 
potential of a sphere with the core density pi and mantle 
density p 2 is 


V(r) = - 


AttGPiRI +P2ir^ -R\) 
3 r 


-2nGp2(Rl-r\ (3) 


where Ri and R 2 are the core and planetary radii (Ri < 
r < R 2 ), respectively, and G is the gravitational con¬ 
stant. The value of (dp/dS)p = -aTp/Cp depends 
on the choices of EOS, V and T but it is typically 
-0.3-0.5 kg^K/m^J over the parameter range of in¬ 

terest. a is the thermal expansion coefficient, and Cp is 
the specific heat at constant pressure. 


2.5. Initial conditions 

Both the target and impactor consist of 70% man¬ 
tle (forsterite or MgSiOs liquid) and 30% core (iron) 
by mass. We assume that the mantle is composition- 
ally uniform prior to the impact (as discussed in Sec¬ 
tion |4^ it is possible that the mantle was composition- 
ally stratified prior to the giant impact, as would be the 


case if there were melt migration and a possible basal 
magma ocean prior to the impact, Labrosse et al.|2007] ). 
Initially, the mantle has a uniform entropy such that 
the temperature at the model surface is approximately 
2000K (3165 J/K/kg for forsterite and 3350 J/K/kg for 
MgSiOs liquid). It should be noted that this parameter 
is not important in the model or in the interpretation of 
the results. Of course, the true physical temperature of 
Earth at its surface could be as low as ~ 300K if any 
earlier steam atmosphere has collapsed. However, the 
high heat flow in this epoch would lead to temperatures 
close to the solidus even at a depth of 10km (far too 
small a distance to be resolved in SPH and a negligible 
contribution to the total heat capacity). The essential 
idea here is that the pre-impact Earth cannot cool effi¬ 
ciently at depth to below the solidus temperature (and 
perhaps not even to the solidus temperature on average) 
because solid-state convection is unable to eliminate the 
heat of previous large impacts by convection in the time 
that has elapsed since the previous giant impact. The 
most important parameter in assessing the melting upon 
impact is irreversible entropy production and not tem¬ 
perature. Here, three types of giant models are investi¬ 
gated: (a) standard, (b) fast-spinning Earth, and (c) sub- 
Earths. The parameters considered are the impactor-to- 
total mass ratio Mi/Mj, the total mass of the target and 
impactor, Mt, the scaled impact parameter b (b = sin 6, 
here 6 is the impact angle), the impact velocity vimp, and 
the initial spin period, Tspin- The initial conditions of the 
models are listed on Table [T] 


3. Results 


3.1 . Comparison of equations of state 

The behaviors of the two EOS are shown in Eigure[2 
One of the differences is that the entropy of the forsterite 
EOS is concave-up, whereas that of the MgSiOs liquid 
EOS is concave down (Eigure[^). This disparity can be 
explained by the signs of {d^S ld^p)T. Eor the MgSiOs 
liquid EOS, this expression becomes {q l)Cvy/p^, 
where y = yo(V/Voy is the Griineisen parameter (dis¬ 
cussed in Section Sl.l) and Cy is the specific heat at a 
constant volume. The exponent q is positive for solids 


and negative for liquids ( [Stixrude et al.||2009] ). There¬ 
fore, the exponent q of the MgSiOs liquid EOS is cho¬ 
sen to be -1.6. On the other hand, the forsterite EOS 
(M-ANEOS) does not differentiate liquids from solids 
and assumes a positive qiy = 7o(WVo)+2/3(l- V/Vq)^, 
Equation 4.11, Thompson and Lauson||1972| ). There¬ 
fore, (d^S/d^p)T is positive in this EOS. The behav¬ 
iors of the internal energies are similar (Eigure[^). The 
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pressures of the two EOS are not significantly different 
at low temperature {T = 4000K) but the difference be¬ 
comes larger at higher temperatures, which stems from 
the difference in q (Figure [^) {{dPldT)p = yCylV is 
proportional to derived from Equation 8, Section 
Sl.l). 


3.2. Mantle structure 

Figure shows a cross section of the Earth’s mantle 
after the Earth reaches its equilibrium state (~ 1-2 
days after the impact). The cross-section lies on the 
z = 0 plane, which includes the center of the planet, 
and is perpendicular to the Earth’s spin axis (z is paral¬ 
lel to the spin axis). The color gradient scales with the 
entropy gain and iron is shown in grey. As discussed 
in Section 2.4 the entropy is higher near the surface 
(dS I dr > 0). Thus, the mantle is thermally stratified (no 
convection). In (a), the impactor hits the surface of the 
Earth twice and becomes disrupted by the impact itself 
and the tide from the planet. The surface of the Earth is 
more shock-heated than the inner part partly because the 
surface is closer to the impact point and partly because 
the heavily shock-heated impactor envelops the surface 
of the Earth. In (b), the impact velocity is so high that 
the impactor penetrates all the way to the CMB (Figure 

bl). This energetic impact strips off outer parts of 
the mantle, which is a desirable feature for producing a 
disk that is composed of the Earth’s mantle rather than 
projectile. This strong impact heats the deep mantle ma¬ 
terials, which are buoyant and come up to the surface, 
whereas the ambient colder materials flow to the inner 
region. Additionally, the outer part of the Earth’s man¬ 
tle is ejected in the z-direction after the impact, which 
eventually falls back and hits the surface of the Earth. 
These materials are highly shock-heated through these 
processes and eventually distributed on the surface of 
the Earth. In (c), the two objects collide each other sev¬ 
eral times and eventually merge into a single planet. The 
mantle becomes highly shock-heated during these pro¬ 
cesses. The surface is more severely shock-heated be¬ 
cause it is close to the impact points. Thus, although 
dS I dr > 0 is a common feature, the reasons for this 
state differ between the models. 

The azimuthal dependence of the entropy is not sig¬ 
nificant. In (a), the impactor hits different parts of the 
mantle. Additionally, the mantle of the impactor, which 
is highly shock-heated, eventually falls on and covers 
the surface of the Earth’s mantle. Through these pro¬ 
cesses, the Earth’s mantle becomes approximately uni¬ 
formly shock-heated. In (b), the mantle shows a slightly 
greater dependence than the other cases. The entropy 
difference at the same r can reach up to a couple of 


hundred J/K/kg, but it has only minor effects on the 
outcome (discussed in Section |4^ . In (c), the mantle 
is approximately uniformly heated through the multiple 
collisions. 

The two EOS provide similar results for each type 
of impact. Figure shows the entropy, density, tem¬ 
perature and pressure of the mantle. We set the mini¬ 
mum density pmin to avoid numerical problems, as de¬ 
scribed in the Section [SL4] In Figure |^, dS jdr > 0 is 
clearly shown. Although the MgSiOs liquid EOS pro¬ 
vides a slightly smaller entropy gain than the forsterite 
EOS, the entropy distribution is very similar between 
the two EOS (Figure [^). A clear difference is that in 
the forsterite EOS, the mantle is “puffier” than in the 
MgSiOs liquid EOS in (b) and (c) (Figure |^, mean¬ 
ing that the forsterite EOS provides a broader region 
with small density (Figure [^). The forsterite EOS has 
a higher temperature at a given r (Figure |^), simply 
because the temperature of the forsterite EOS needs to 
be higher to explain the same entropy as the perovskite 
EOS in the density range examined (Figure [^). The 
pressure distributions are very close until reaching the 
outermost part (Figure |^). Nevertheless, the overall 
physics of the two mantles are not significantly differ¬ 
ent. In addition to this, we also perform a simulation 
for (a) with an MgSiOs bridgmanite EOS and show that 
this result is similar to that of the MgSiOs liquid EOS 
(Section [S1.3| ). 

3.3. Mantle melting and mixing 

Figurej^ shows that the majority of the mantle expe¬ 
riences melting (A^meit = 623 J/K/kg). Specifically, the 
fractions of the mantle that are melted are (a) 68%, (b) 
99 %, and (c) 100% for the MgSiOs liquid EOS and (a) 
80%, (b)100%, and (c) 100% for the forsterite EOS (if 
A^meit = 500 J/K/kg, the fractions become (a) 80%, (b) 
100%, and (c) 100% for the MgSiOs liquid EOS and (a) 
92%, (b)100%, and (c) 100% for the forsterite EOS. If 
A^nieit = 200 J/K/kg, the fractions become (a) 98% for 
the MgSiOs liquid EOS and the other cases are 100%). 
Therefore, even if there is no magma ocean prior to the 
impact, the base of the magma ocean is close to CMB in 
all cases. The analysis of the mantle mixing is summa¬ 
rized in Tablej^ The potential energy is normalized by a 
constant (0.01 x where Vesc is the escape veloc¬ 

ity of the Earth). A^ave represents the average entropy 
gain due to the impact. In (a), the ratio APE/AKE is 
much greater than 0.5 in both EOS, which indicates that 
the mantle may remain unmixed. In (b), the ratio is still 
above the critical value, but not as significantly. Given 
the uncertainties in our model, it is difficult to deter¬ 
mine the stability of the mantle in this case. This would 
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presumably indicate that the mantle or at least part of 
the mantle may remain unmixed. An important factor 
here, however, is that once an additional kinetic energy 
available to mix the mantle is considered, it becomes 
more likely that the mantle becomes mixed in (b) (dis¬ 
cussed in Section |43| ). In (c), the ratio is much less than 
0.5, which may lead to dynamical mixing of the man¬ 
tle. Since AKE has larger variations than APE among 
the models, the determining factor is the kinetic energy. 
The kinetic energies of the impacts determined by the 
different models rank as follows: (a) < (b) < (c) (Table 
[^. It should be noted that although the endpoints of the 
kinetic energy of the mantle in (b) and (c) are similar, 
AKE in (b) is significantly lower than that in (c) because 
the pre-impact Earth has a large amount of kinetic en¬ 
ergy due to the rotation. Thus, our study indicates that 
the standard model (a) is most consistent, the sub-Earths 
model (c) may be least consistent with the preservation 
of the mantle heterogeneity, and the fast-spinning Earth 
(b) may lie between. 


4. Discussion 


4.1. A discrepancy between SPH and the Rankine- 
Hugoniot equations 

The Rankine-Hugoniot equations describe the rela¬ 
tions between pre- and post-shock states of materials 
(e.g., |Tonks and Melosh |1993t [Stewart et al.[ |2014| ). 
The equations can predict the entropy increase of a 
shocked material by a small impact, but these equations 
do not model a giant impact very well. This is because 
the entropy increases by a giant impact not only be¬ 
cause of the primary impact-induced shock, which pro¬ 
vides the peak shock pressure (as shown in Eiguresj^al- 
cl and 
Eigure 


S2] ), but also because of subsequent processes, 
^hows a discrepancy between the Rankine- 
Hugoniot equations and SPH calculations. It exhibits 
the entropy of each SPH particle (shown in a dot) af¬ 
ter the impact as a function of the shock peak pressure. 
The lines represent predictions of entropy gain using 
the Rankine-Hugoniot equations ( [Sugita et al.|2Q12 


see 


Section |S1.5| ). It is shown that the entropy gain pre¬ 
dicted by SPH is generally much larger than the predic¬ 
tions by the equations, especially in (b) and (c). This is 
because the mantle is heated due to a number of addi¬ 
tional processes, such as reflected shock waves within 
the Earth as well as planetary deformation followed by 
gravitational potential energy release. The latter effect 
is clearly shown in Eigure which represents the pres¬ 
sure evolution during the giant impact. In (b) and (c), 
after the impact, a large fraction of the Earth’s mantle 


experiences an extensive expansion and becomes sub¬ 
ject to low pressure (b2 and c2). Subsequently, the part 
falls back towards the Earth’s core and gains a high pres¬ 
sure (b3 and c3). Through these processes, the part of 
the mantle efficiently gains entropy. In (a), the extent of 
such deformation is not as prominent as the other cases; 
therefore, its entropy increase is closer to the analyti¬ 
cal estimates. The time dependent relationship between 
the pressure and entropy is discussed further in Section 
|S1.6[ Thus, the Rankine-Hugoniot equations of state 
predict only a part of the entropy gain (due to the pri¬ 
mary impact-induced shock); therefore they do not rep¬ 
resent the total entropy gain by a giant impact. 

4.2. Mixing analysis 

Eor this mixing analysis, we assume that AKE is the 
kinetic energy of the differential rotation of the mantle. 
However, differential rotation is not the only source of 
energy available for mixing. In addition to flows pro¬ 
duced by the impact itself, there will be excitation of 
normal modes of the planet (we have observed these in 
longer runs of SPH simulations for close encounters or 
glancing impacts). These are analogous to waves on an 
ocean and have an energy that oscillates between grav¬ 
itational energy and kinetic energy of the fluid motion, 
with equipartition of the average kinetic energy and av¬ 
erage gravitational energy. The mean kinetic energy is 
therefore simply related to the gravitational energy for 
a distortion of an equipotential surface or the planet as 
a whole by a vertical distance h. These normal modes 
will “ring down” slowly, but in a region of static stability 
the distortion of equipotential surfaces creates horizon¬ 
tal density gradients that drive rapid baroclinic instabil¬ 
ities. Thus, some significant fraction of this energy (the 
precise amount is not easily estimated) can go into mix¬ 
ing. Here we make a simple estimate for the additional 
kinetic energy for mixing: the gravitational energy by 
distortion for a uniform density sphere is proportional to 
^^(^)^. Approximately, is ~ 0.1 in (a) and > 0.3 
in (b) and (c). Thus, the gravitational energy becomes 
~ 1, 10 and 10 (in 0.01 which is the unit of 

AKE in Table 1^. Assuming that approximately half of 
this energy eventually goes into heat and the other half is 
available to mix the core (this increases AKE), corrected 
APE/AKE become 1.60, 0.35, and 0.17 in (a), (b) and 
(c) (MgSiOs liquid EOS). Thus, the mantle would still 
remain unmixed in (a), but the mantle would be mixed 
in (b) and (c). 

We ignore the azimuthal dependence, which is likely 
a reasonable approximation. To demonstrate that the 
approximation is feasible, consider the worst case, the 
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mantle in (b), whose entropy distribution shows the 
strongest azimuthal dependence. We evenly divide the 
mantle into four sections according to its azimuthal an¬ 
gle. We observe that the average APE in each section 
differs from the globally averaged value by less than 
10%, which is not sufficient to alter the mixing status. 
Another potential concern is that we only perform one 
simulation for each model and EOS, therefore it is pos¬ 
sible that the ratio APE/AKE varies among models, but 
we show that the ratio we derived is likely to capture the 
general trend (discussed in Section |SE7]). 


It should be noted that the mixing analysis described 
in Section only determines the global mixing state. 
Therefore, it may be possible to preserve a local het¬ 
erogeneity even at APE/AKE < 0.5, especially if the 
ratio is close to the critical value. Unfortunately, our 
SPH does not provide information about local mixing. 
To investigate such local mixing, it would be necessary 
to perform a simulation with much higher resolution, 
but doing so would be computationally quite expensive. 
Additionally, we do not consider the possibility that the 
mantle was compositionally stratified before the impact 
(e.g. a denser layer at the base of the mantle), but if 
this was the case, it could increase the stability of the 
mantle. Moreover, we discuss issues and implications 


of this mixing analysis in Section SI.7 


4.3. Mixing during subsequent cooling 

We turn now to the question of whether the mantle 
will mix during the post-giant impact cooling. This is 
a simple problem, conceptually, and it does not require 
extensive numerical analysis. As the outer regions of 
Earth cool, the mantle will evolve into two regions, a 
nearly isentropic region that can convect and a deeper 
region where the entropy is still lower than that appro¬ 
priate for the effective radiating temperature of Earth at 
that time. The deeper region that is of lower entropy 
is unable to convect since it is still part of the stably 
stratified zone. Thus, the convective zone propagates 
downwards as the planet cools and the deepest part of 
the mantle is only able to convect when the entire re¬ 
gion above has cooled to the entropy that this deepest 
part had immediately after the giant impact. But as we 
have seen, this entropy is not much larger than it was 
before the giant impact, at least in the case of the canon¬ 
ical giant impact. This takes a long time relative to the 
time it takes to make the Moon, as we now demonstrate. 
Let Te be Earth’s effective temperature for the thermal 
state corresponding to the isentrope that is only slightly 
higher (and hotter) than the pre-giant impact state. At 


that epoch, the cooling equation takes the form 


= ~{MCvT^) 

dt 


( 4 ) 


where M is the mantle mass and is the mean temper¬ 
ature at that epoch. Accordingly, the time r that takes to 
cool an amount Ar^ is 


T ~ 5 X 10' 




1000K/\ T, 


pOOOKV 


yrs 


(5) 


This is an underestimate because we assume a typical 
specific heat rather than the additional term that comes 
from freezing. Importantly, Te can be low at this point, 
though its exact value depends on whether there is a 
steam atmosphere. Eor example, at Te= 500K, this cool¬ 
ing time for lOOOK is tens to hundreds of thousands of 
years, which is much longer than the estimated time for 
making the moon (hundreds to thousands of years). The 
reason for this is obvious: Earth’s mantle has a much 
greater heat capacity than the Moon-forming disk and 
yet a smaller area from which to cool. This is the time 
that must elapse before any mixing of the deepest man¬ 
tle material can even begin. Even at that time, mixing 
will not occur if there is a basal magma ocean prior to 
the impact. It is often assumed in the literature that mix¬ 
ing occurs if there is liquid but it is important to stress 
that the state of the material is not the crucial issue here. 
Counterintuitive though it may seem, it is harder to mix 
a stably stratified liquid than to mix a solid of the same 
stratification because the thermal anomalies associated 
with convection in a liquid are much smaller than those 
for solid state convection. 


4.4. Impactor's isotopic signature 

Eor the standard scenario (a), the impactor needs to 
have isotopic signatures similar to those of the Earth 
to explain the geochemical evidence. To gain a better 
idea of how likely this similarity could be in the Grand 


Tack model, we run similar analyses to those of Pahle- 


van and Stevenson ( 2Q07| ). Although |Kaib and Cowan 
( |2015 1 have conducted isotopic analyses of a model that 
is similar to the Grand Tack, the initial conditions are 
more simplified than the original model. 

We take 30 planets from eight simulations from 
Walsh et ah] ( |2011| ) (each simulation produces three to 


four terrestrial planets) and examine the difference in 
the oxygen isotope ratios between the target and im¬ 
pactor. We assume that the oxygen isotopic ratio in the 
protoplanetary disk linearly varies as a function of the 
heliocentric distance r', A^^O(r') = ci x r' -l- C 2 . Two 
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pairs of ci and C 2 are used; case 1: ci and C 2 are chosen 
such that Ai^O(r') = 0 %o at 1 AU and A^^O(r') = 0.321 
%c at 1.5 AU(A^^O of Mars), and case 2: ci and C 2 are 
chosen such that the Earth and Mars analogues in each 
simulation have A^^O(r') = 0 %c and 0.321 %c, respec¬ 
tively (this is the same setting as those of |Pahlevan and 


|Stevenson|2007| ). Simulations that do not produce Mars 
analogues are discarded. Figure shows histograms of 
the differences in the oxygen isotopic ratios. The top 
and bottom planes correspond to case 1 and case 2. The 
left panels include all of the impacts and the right pan¬ 
els focus on the last giant impact for each planet. Even 
if we take the stricter criterion (the isotopic difference 
should be less than 0.005 %c), there is one last giant im¬ 
pact out of 30 events and one out of 19 events that satisfy 
this criterion in case 1 and case 2, respectively (Figure 
[fright panel). It should be noted that the mass involved 
in these particular impacts are a couple of Mars masses, 
so that these do not mimic the Moon forming impact. 
Nevertheless, this simple analysis may indicate that the 
Grand Tack model could produce an impactor that has 
a similar composition to the Earth and its possibility is 
around 3-5 %. Needless to say, a much larger dataset is 
needed to estimate the probability. Additionally, other 
isotopic ratios, such as tungsten and silicon, need to be 
explained in addition to oxygen (Section [^. 


5. Conclusion 

To investigate the relationship between the Moon¬ 
forming model and geochemical evidence, we have in¬ 
vestigated the initial state of the Earth’s mantle after the 
Moon-forming impact. Giant impact simulations using 
the SPH method are performed with the MgSiOs liq¬ 
uid and forsterite equations of state. Three impact mod¬ 
els are considered here: (a) standard, (b) fast-spinning 
Earth, and (c) sub-Earths. We find that the mantle be¬ 
comes mostly molten after the Moon-forming impact 
in all of the cases. It is also shown that the impact- 
induced entropy gain of the mantle cannot be predicted 
well by the Rankine-Hugoniot equations. Based on our 
analysis on the energy balance of the mantle, we find 
that the Earth’s mantle is likely to stay unmixed in (a), 
while it may be mixed in (b) and (c). The extent of 
mixing is most extensive in (c). This is primarily be¬ 
cause the impact-induced kinetic energy in (c) is much 
larger than that of (a) while (b) falls between. Thus, the 
standard model (a) is most consistent, the sub-Earths 
model (c) may be least consistent with the preserva¬ 
tion of the mantle heterogeneity, and the fast-spinning 
Earth (b) may lie between. We also find that if the im¬ 
pact does not mix the mantle, the Earth’s mantle would 


remain partially unmixed for more than the Moon ac¬ 
cretion time scale (~ 100 - 1000 yrs). It is therefore 
possible, at least in the standard giant impact case (and 
possibly some others), to avoid complete mixing during 
and immediately after a giant impact. The issue of sub¬ 
sequent mixing (>~ 1000 yrs) is not a simple problem 
but is no different from the usual assessment of mantle 
convection mixing (e.g., Tackley| 2012| ). 
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Figure 1: The two equations of state at various temperatures. The blue, orange and red lines correspond to 4000K, 8000K, and 15000K, respectively. 
The solid and dashed lines represent the MgSiOs liquid and forsterite EOS. 
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Forsterite MgSiOs liquid 



Figure 2: Cross-section of the mantle after the impact. The color gradient scales with the entropy gain of the mantle material. Iron is shown in 
grey. The top panels show the case with the MgSiOs liquid EOS and the bottom panels show the case with the forsterite EOS. 
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Figure 3: Mantle structure as a function of r. The blue, orange and grey lines represent the standard, fast-spinning and sub-Earth cases. The 
solid and dashed lines represent the MgSiOs and forsterite EOS. The dash-dot line represents the entropy of melting. Although the density and 
temperature are slightly different, the entropy is quite similar for the different EOS. The pressure is also nearly the same except for the outer 
boundary. 
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Figure 4: Entropy gain of each SPH particle in the mantle with the MgSiOs liquid EOS. The dashed (pi = 0 GPa) and solid (pt = 50 GPa) lines 
represent estimated entropy gain based on the Rankine-Hugoniot equations, where pt is the pre-shock pressure. To examine the lower part of the 
mantle, this figure only shows SPH particles that were initially under high pressure (> 24 GPa) before the impact. 
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Figure 5: Snapshots of early stages of the impacts (al-3: 0.62, 1.0, 1.5 hrs, bl-3: 0.37, 0.60, 0.70 hrs, and cl-3: 0.33, 0.60, 0.72 hrs after the start 
of the simulation). The color gradient scales with the pressure of the mantle material and iron is shown in grey. The Earth deforms significantly 
and experiences low pressure in (b) and (c) while the deformation is less intense in (a). 
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Figure 6: Histogram of the differences in the oxygen isotopic ratios between the target (proto-Earth) and impactor. The left panel includes all of the 
giant impacts and the right panel includes only the last giant impact for a planet. The top and bottom panels represent different initial conditions 
(case 1 and case 2, Sectionj^J. The dashed and dash-dot lines indicate the two estimates of the oxygen isotopic differences between the Earth and 
Moon (0.005%o and 0.012%o, respectively, |Wiechert et al.|2001||Herwartz et al.|2014) . 
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MJMt: 

Mj 

b 

rimp 

^spin 

(a) Standard 

0.13 

1.02 

0.75 

1.0 

0 

(b) Fast-spinning Earth 

0.045 

1.05 

-0.3 

20 (km/s) 

2.3 

(c) Sub-Earths 

0.45 

1.10 

0.55 

1.17 

0 


Table 1: The initial conditions of the impact. M^IMj is the impactor-to-total mass ratio, Mj is the total mass scaled by the mass of the Earth, b 
is the scaled impact parameter, vimp is the impact velocity scaled by the escape velocity (except b), and Tspm is the initial spin period of the target 
(hrs). 



AKE 

APE 

APE/AKE 

^S ave 

(al) Standard, liq 

1.38 

3.01 

2.18 

1181 

(a2) Standard, for 

1.51 

2.82 

1.86 

1220 

(bl) East-spinning Earth, liq 

3.08 

2.79 

0.91 

1645 

(b2) East-spinning Earth, for 

3.36 

2.45 

0.73 

1675 

(cl) Sub-Earths, liq 

5.85 

1.85 

0.32 

2210 

(c2) Sub-Earths, for 

5.50 

1.26 

0.23 

2196 


Table 2: The outcomes of the impact. AKE is the change in kinetic energy (normalized by 0.01 x where Vesc is the escape velocity of 

the Earth) and APE is the potential energy required to mix the mantle. A*Save is the average entropy gain of the mantle, liq represents the MgSiOs 
liquid and for represents forsterite EOS. 


SI. Supplementary Materials 


The entropy 5 (V, T)is described as 


Sl.l. Equation of state 

Our equation of states for the MgSiOs liquid fol¬ 
low the formulae previously derived (e.g., |de Koker and 
|Stixrude[ |2009[ [Stixrude et al.[ |20Q9| ). The Helmholtz 
free energy is written as 

F{V, T) = Fo + Fcmp(V, To) + Fth(V, T), (1) 

where Fq = F(Vo, To) is the free energy at the refer¬ 
ence volume, Vo, and temperature. To. T’cpmCV, To) and 
Fth(K T) are the compressional and thermal contribu¬ 
tions to the free energy, respectively. Fcmp is 


Fcmp=9^roVo|^/' + ^a3/|, (2) 

where 

03 =3(/i:^o-4), (3) 

f=\ - ij. (4) 

Kto is the isothermal bulk modulus (at T = Tq), and 


K'r^Q is its pressure derivative at /? = 0 and T = To. T’th 
is written as 


/^th 



s{yx)dr. 


(5) 


S(V,T) = S 


JVo 


Cv7(V', To) 


V 


dV' 


f 

OTo 


Cy(v, r) 


dr, 


( 6 ) 


where So = S(Vo,To), Cy is the specific heat that is 
assumed to be a constant, y is the Griineisen parameter, 
which is described as 


r = ro 



(7) 


where, yo and q are constants. 

The internal energy, E, and pressure, P, are described 
as 


E{V, T) = Eo + 9 KtoVo I If + Usf 


-f Cv(T - To) + CyTo 


f r(v\ 


To) 


dV, 

( 8 ) 


P(y, T) = -iKToH + + ^f) 

y(v, To) 

^Cy(T-To)^-^^. (9) 

Here, Eo = E{Vo,To). po = (^IVo)^ To, Kjo, K^o^ Cy, 
yo, q, Eo, and are listed in Table 
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Po(kg/m3) 

To(K) 

Xro(GPa) 

KjQ 

Cv(J/K/kg) 

To 

4 

£o(MJ/kg) 

5o(kJ/K/kg) 

2650 

2000 

27.3 

5.71 

1480 

0.6 

-1.6 

2.64 

3.33 


Table SI: Constants for the MgSiOs liquid EOS. 


SI.2. Mixing criterion 

In a simple shear flow, the criterion for a Kelvin- 
Helmohltz instability (effectively the criterion for mix¬ 
ing) is Ri= N^liduldz)^ < 1/4. Ri is the Richardson 
number of the system, N is the Brunt-Vaisala frequency, 
= -g{dpldz)lp, u is the velocity, z is the direction 
perpendicular to the flow, g is the gravity, and p is the 
density. Ri is related to the ratio of potential energy 
to kinetic energy (e.g., Taylor 19311 Chandrasekhar 


1961) and is normally defined in terms of fluids with 
well defined constant density differences or a density 
gradient, but our system has variable values of all the in¬ 
put parameters. Therefore, we must necessarily restate 
the problem in terms of the energy budgets rather than 
explicitly in terms of velocity shear. The kinetic energy 
difference AKE per unit area for a layer of thickness L 
between the initial state (with shear) and the final state 
(with no shear but the same linear momentum) is 

1 

-p (u{zf-iuol2f)dz = Lpull24, (10) 

^ J-L/2 

where uq is the initial velocity difference between the 
top and bottom and p is the mean density. The gravita¬ 
tional potential energy difference (APE) between initial 
and final (fully mixed) states is 




r 

I (P - p(z))gzdz = pN^C /12, 

J-LI2 


( 11 ) 


where the density variation is assumed small. There¬ 
fore, the Richardson number criterion becomes 
AKE>2APE and the physical interpretation is that one 
must provide not only the energy to overcome the poten¬ 
tial energy difference but also the energy to mix (which 
shows up as heat from the dissipation of small scale 
turbulent motions). In the analysis provided by [Chan 
drasekhar ( 1961| ) (p. 491), the APE he defines is for 
complete overturn (that is, the new density profile is the 
exact opposite of the initial density profile), which is a 
different setting from ours. Therefore, his criterion dif¬ 
fers from ours (mixed if AKE > APE in his anaysis). 

SI.3. MgSiO^ bridgmanite EOS 

Eigure shows cross-sections of the mantles after 
the impact with the MgSiOs liquid and MgSiOs bridg¬ 
manite EOS. The thermodynamic parameters for the 


bridgmanite EOS is listed in Table The entropy 
gains are slightly different, but the extent of shock¬ 
heating and the feature of dS I dr > 0 are similar among 
these cases. One might expect that a liquid mantle may 
gain higher entropy than a solid one based on the study 
done by |Karato] ( |2014| ) that suggests that the surface of 
a molten mantle gains higher entropy by impact than 
a solid surface due to its smaller sound speed Co and 
negative q for the liquid (for the definitions of these pa¬ 
rameters, see Section [ST3] ). However, the difference of 
the two EOS in Co becomes smaller at a greater depth. 
This would diminish the difference in the entropy gain 
among the two EOS. In addition, Karato ( 2014| ) use the 
Rankine-Hugoniot equations to describe the physics of 
the planetary surface, but we cannot use these equations 
to predict the entropy gain of the entire mantle as dis¬ 
cussed in Section 4.1. 


SI. 4. Correction of the outer boundary 

The density of the outer edge of the mantle is cor¬ 
rected because the simulation itself does not provide an 
accurate value. One of the reasons is that the standard 
SPH cannot describe a large density difference (e.g., 
CMB or planet-space boundary). The density of a par¬ 
ticle at the outer boundary becomes too small because 
the particle does not have many nearby particles; thus, 
the smoothing length becomes large. This leads to a 
problem that dPjdr at the outermost part of the mantle 
becomes nearly 0 or it even becomes positive in (b) and 
(c) (Eigure 3D). This state is not physically sensible be¬ 
cause the hydrostatic equation is not correctly solved in 
the region. 

To avoid this numerical problem, we define the min¬ 
imum density pmin = p(^max)- Here, r^ax is the maxi¬ 
mum r whose region satisfies dPjdr < 0. If the den¬ 
sity at r > Tmax is lower, the r is recalculated by set¬ 
ting p = pmin and conserving the mass. Typically, 
Pmin ~ 1500 - 1600kg/m^ (Eigure 3B). This is uncer¬ 
tain but may be reasonable because this is close to a 
rough estimate of the density at such a high tempera¬ 
ture. The density at the outer edge can be approximated 
asp ~ po(l- aT) ~ 1577kg/m^ ata = 2.7x10“^ ( Eiquet 


et al. 


2000| i, po = 2650kg/m3, and 7 = 1.5 x lO^K. Af- 


ter this procedure, APE in the two EOS become similar. 
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Figure SI: Entropy of the Earth’s mantle after the impact. The top panel shows the case with the MgSiOs liquid EOS (the same as Figure 2 in the 
main text) and the bottom panel shows the case with the MgSiOs bridgmanite EOS. 
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Po(kg/m^) 

To(K) 

/fro(GPa) 

KjQ 

Cv(J/K/kg) 

7o q £o(MJ/kg) 

5o(kJ/K/kg) 

3680 

2000 

200 

4.14 

1200 

1.0 1.0 1.995 

2.63 


Table S2: Constants for the equation of state of MgSiOs bridgmanite. 


Thus, although this approximation is simple, it provides 
a reasonable answer. 


SI.5. The Rankine-Hugoniot equations 


[Sugita et~n? ( 2012 1 derive the following differential 
equations to describe after-shock temperature T and en¬ 
tropy S based on the Rankine-Hugoniot equations; 


dT AUs-U.T 

— = CoroT- 

dUp u 


q+l 


sUl 

CvUs 


( 12 ) 


dS _ sUl 

^ " tW; 


(13) 


Here, p — pi pfUgUp and p — piJJgjiJJg JJp\ 
where pi, pi. Us, and Up are the pre-shock pressure, pre¬ 
shock density, shock velocity and particle velocity. Us 
and Up has a relation Us = Cq + sUp, where Cq and 5- 
are the sound speed and constant. For our calculations, 
we choose 5' = 1.56 (for MgSiOs bridgmanite, |Deng| 
et al.||2008] ) and Tt = 2000K (pre-shock temperature). 
At Pi = 0 GPa, pi = 4100kg/m^, Co = 6.47km/s and at 
Pi = 50 GPa,p, = 4500kg/m^ Co = 9.0km/s. 


SI.6. Pressure entropy increase 

Figure shows the relationship between the pres¬ 
sure (shown in grey) and entropy gain (shown in green). 
We choose a specific SPH particle from each simula¬ 
tion and follow its properties. In (a), the primary impact 
(~ 90 GPa) is the major source for the entropy increase. 
The entropy changes overtime, but the extent is limited. 
In (b), the SPH particle is heated by multiple shocks, in¬ 
cluding the primary impact-induced shock and shocks 
due to the planetary expansions and contractions (dis¬ 
cussed in Section 4.1). After ~ 5 tirs, the entropy slowly 
increases due to continuous small-scale planetary defor¬ 
mation (the planet continues to wobble) until the system 
reaches its equilibrium state. In (c), the SPH particle 
experiences a number of shocks because the target and 
impact collide several times. The entropy gain is larger 
than the other two cases. 


SI. 7. Further discussions on the mixing analysis 


We assume that the Earth’s mantle was chemically 
heterogeneous before the impact, but here we further 
discuss its plausibility. Unlike (a) or (c), the model (b) 
requires that the Earth spins very quickly before the gi¬ 
ant impact. This may indicate that the Earth experienced 
another (older) giant impact before the Moon-forming 
impact. This is because the angular momentum of a 
planet delivered by a number of small impacts from ran¬ 
dom directions tend to cancel out. This older giant im¬ 
pact could have been similar to the “sub-Earths” model, 
meaning that two similar mass objects collided, because 
this type of an impact is one of the easiest ways to de¬ 
liver a large angular momentum to the planet ( |Canup| 
|2014| ). If this is the case. Earth’s mantle could have 
been homogenized before the Moon-forming impact. 
If the heterogeneity formation predated this older im¬ 
pact, this could be a potential problem for (b). Alterna¬ 
tively, it is also possible that the heterogeneity formed 
between this older and the last giant impacts, possibly 
in the form of a basal magma ocean by fractional melt¬ 
ing and crystallization processes. The re-establishment 
of a compositionally distinct basal magma ocean could 
have been accomplished in less than 10^ yr compared to 
the likely time between giant impacts plausibly ~ 10^ 
years (without an atmosphere, the majority of the man- 
tle could have crystallized as short as 10^ years, [Sola 


|matov||2000| ). It should be noted that there is no rea¬ 
son to suppose that this older impact was immediately 
prior, but the time interval might have been preferably 
short if the older impact formed a satellite. This is 
because the interaction between the Earth and satellite 
may have slowed the Earth’s spin rate within 10^ - 10^ 
years ( jCanup 2014 ). This older satellite might have 
merged with a newer satellite formed by the last giant 
impact ( [Citron et^|2014| ). 

Another potential problem is that the critical value 0.5 
has been derived to analyze the stability of horizontally 
stratified layers, but the value can differ for spherically 
stratified layers, as in our model. However, there is no 
literature concerning this geometry; thus, we simply ap¬ 
ply the critical value for our analyses. The choice of the 
minimum density could also affect the estimate of APE. 

Furthermore, we only perform one simulation for 
each model and EOS. It is possible that AKE and APE 
can change even for the same type of impact depend- 
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Figure S2: Time dependence of the pressure and entropy (in 10^ J/K/kg) of an SPH particle during the initial 24 hours. 


20 













ing on the choice of the initial conditions (e.g., Vimp 


and 


b). To perform a simple and crude analysis, here we as¬ 
sume that the planetary kinetic energy is expressed as 
where / is the moment of inertia and oj is the 
angular velocity of the planet, and that 1 and APE do 
not vary in the same model. We compute the ratio of 
APE/A KE based on oj from pubhshed success ful simu¬ 


lations {Cuk and Stewart 


2012 


Canup 


2012) and find 


that most of these simulations do not change the ra¬ 
tio large enough to cross the critical value 0.5, except 
Run 14 (Mi/M t= 0.45, ^7=0.40, and Vimp/Vesc = 1-4) in 
the sub-Earths model, depending on the EOS (0.52 for 
MgSiOs liquid and 0.38 for forsterite). Thus, our re¬ 
sults likely provide the general trend, but some outlier 
may exist. Nevertheless, the choice of initial conditions 
is not likely to alter the signatures of dS I dr > 0 or the 
melting of the nearly entire mantle because these are 
states less sensitive to the conditions. 
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